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ABSTRACT 

In X-ray binaries, several percent of the compact object luminosity is intercepted by the surface of the normal 
companion and re-radiated through Compton reflection and the K-fluorescence. This reflected emission follows 
the variability of the compact object with a delay approximately equal to the orbital radius divided by the speed 
of light. This provides the possibility of measuring the orbital radius and thus substantially refining the compact 
object mass determination compared to using optical data alone. We demonstrate that it may be feasible to measure 
the time delay between the direct and reflected emission using cross-correlation of the light curves observed near 
the Ka line and above the K-edge of neutral iron. In the case of Cyg X-1, the time delay measurement is feasible 
with a 3 X 10^- 10'' s observation by a telescope with a 1000 cm^ effective area near 6.4 keV and with a ~ 5eV 
energy resolution. With longer exposures, it may be possible to obtain mass constraints even if an X-ray source in 
the binary system lacks an optical counterpart. 

Subject headings: binaries: spectroscopic — X-rays: general — ^X-rays: stars 



L INTRODUCTION 

The existence of stellar-mass black holes in our Galaxy is al- 
most exclusively established by measuring the mass of the com- 
pact object in several bright X-ray binaries. A compact object 
more massive than approximately SM© is unstable and should 
collapse to form a black hole (Rhoades & Ruffini 1974, Chitre 
& Hartle 1976). Therefore, if the dynamically measured mass 
of a luminous X-ray source exceeds 3Mq, it is considered a 
firm black hole candidate. 

In X-ray binaries, the orbital velocity of the compact object 
cannot be measured by means of optical spectroscopy. This 
means that the only directly measured quantity is the mass func- 
tion /(M), 



/(M) = 



Mi sin / 



(1) 



(M,+Mopt)2' 

where M^ and Mopt are masses of the compact object and normal 
companion, respectively, and / is the inclination angle of the or- 
bit. The mass function is the lower limit of the compact object 
mass. In many high mass X-ray binaries, the mass function is 
uninterestingly low; for example, f{M) = 0.25 Mq for Cyg X-1 
(Gies & Bolton 1986). A determination of Mv is possible if Mopt 
and sin/ are independently constrained. Extensive work has 
been done to estimate these parameters for several black hole 
candidates (Bolton 1972, Liutyj, Sunyaev & Cherepashchuk 
1975, McCHntock & Remiflai-d 1986, Remillard, McClintock 
& Bailyn 1992, HeiTero et al. 1995, Filippenko, Matheson & 
Barth 1995, Bailyn et al. 1995, Remiflard et al. 1996). How- 
ever, the derived compact object mass can be rather uncertain 
due to its strong dependence on Mopt and especially sin/. For 
Cyg X-1, mass values in the range 6-I6M0 are found in the 
literature (Dolan & Tapia 1989, Gies & Bolton 1986). 

The mass measurement will be much more accurate if some 
information about the orbit of the compact object is available. 
Unfortunately, very fast motions of the material in the vicin- 
ity of the black hole make it impossible to measure the orbital 
velocity from the Doppler shifts of spectral lines. Instead, we 
propose to measure the orbital radius using the following ap- 
proach. In high mass X-ray binaries, some fraction of the com- 



pact object X-ray luminosity is intercepted by the normal com- 
panion surface (the intercepted fraction in low mass binaries 
is too low for an application of our method). A significant 
fraction of this energy is re-radiated in the X-ray band either 
through Compton reflection or line fluorescence. This emission 
reaches the observer with a time delay roughly equal to a/c, 
where a is the separation between the normal companion and 
the compact object. Since the compact object emission is usu- 
ally variable, it may be possible to measure the time delay and 
thus estimate the orbital radius. We show below that this mea- 
surement can be carried out using temporally resolved high res- 
olution spectroscopy of the fluorescent Ka line of neutral iron. 
Let us demonstrate how this new information refines the mea- 
surement of the compact object mass. Suppose for simplicity 
that the orbit is circular and the radius of the normal compan- 
ion is negligible compared to the orbital radius. The time delay 
at orbital phase (j) is 

Af = fl/c[l-sin(/cos27r(?!))]. (2) 

Equations ([l|)-(g) together with the Kepler's third law, 

fl3 = p2^M, +Mopt)G/47r^ (3) 

can be used to express M^ as a function of just one parameter, 
Mopt, ', or a. In particular, if the time delay measurement is 
performed at inferior conjunction (0 = 0.5), 



M,: 



sin / ( 1 + sin i)^ 



or Ml : 



(M,+Mopt)/'/^ 

gV3_(M, + Mopt)V3^ 



(4) 



where g = 4Tr^(cAt)^/P^G. If the time delay is measured at 
several orbital phases, all system parameters can be determined 
independently. Moreover, useful mass constraints are possible 
in this case even if no optical data are available (§ Q). 

2. MEASUREMENT STRATEGY 

The X-ray emission reflected from the surface of the nor- 
mal companion consists of the Compton-reflected continuum 
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and the fluorescent emission. While the reflected continuum is 
weak compared to the direct emission from the black hole, the 
fluorescent line emission can be relatively strong. The strongest 
fluorescent line is the iron Ka doublet (Basko, Sunyaev & 
Titarchuk 1974). Basko (1978) estimates the equivalent width 
of w 7eV for this line in Cyg X-1. With « 5eV spectral reso- 
lution of X-ray calorimeters, it should be possible to measure 
the source light curve in a narrow spectral interval around Kai 
and K0C2, where the fluorescent line contributes w 50% of the 
total flux. The time delay can be measured by cross-correlating 
the light curve in the Ka region with the continuum light curve 
which is dominated by the direct black hole emission. 

Ka lines are primarily excited through photoelectric absorp- 
tion in the K-edge (Eg = 7.1 keV for neutral iron). There- 
fore, the energy band just above 7.1 keV is best for the con- 
tinuum light curve measurement. The exact choice of the en- 
ergy boundaries is unimportant as long as the source variability 
patterns do not change with energy, as is the case in Cyg X-1 
(Nowak et al. 1999). Below, we assume that the continuum is 
measured in the 7. 1-9.1 keV energy band. 

The most challenging part of the time delay measurement 
is determining the light curve in the Ka line region. This 
Une consists of two components, Kai and Kaa, with energies 
Eac, = 6.404 and £„, =6.391 keV (Bambynek et al. 1972). Both 
components have approximately equal natural widths of 3.5 eV 
(FWHM). The Doppler broadening of the lines due to thermal 
motions of the emitting atoms AE^h = 0.1(7/ 10^ K)'/^ eV is 
much smaller than the natural width of the lines. Therefore, 
more than 90% of the total Ka line flux is within a narrow, 
20 eV, energy interval around E = 6.398 keV. Selecting pho- 
tons in this energy interval is easy with the 5 eV (FWHM) en- 
ergy resolution of future X-ray calorimeters. An approximately 
20 eV energy interval seems optimal for the time delay mea- 
surement. A wider band would include a larger contribution 
of the direct continuum without increasing the line flux signifi- 
cantly. Narrower intervals would miss flux in the line wings. 

The cross-correlation function of the continuum and line 
light curves should contain two prominent peaks. A strong 
narrow peak around the delay r = corresponds to the auto- 
correlation of the emission coming directly from the black hole. 
Reflection from the normal companion results in a secondary 
broad peak in the correlation function. This peak should have a 
centroid at r w -a/c and width At w R/c, where a and R are 
the radii of the orbit and normal companion, respectively. The 
amplitude of the secondary peak should be proportional to the 
equivalent width of the reflected Ka line and the rms^ of the 
source flux variability. Therefore, the time delay measurement 
is best performed when the reflected line is strong (i.e. near in- 
ferior conjunction, Basko 1978), and when the source is highly 
variable. 

In the next section, we simulate an observation of the time 
delay in Cyg X- 1 which could be performed by a telescope with 
an effective area of 1000 cm^ in the 6-9 keV energy range. Such 
an instrument approximately corresponds to one spacecraft of 
the proposed Constellation-X mission.' 

3. SIMULATED OBSERVATION 

We adapt the binary system parameters from Herrero et al. 
(1995) — a radius of the normal companion R = 11 Rq, a sep- 
aration between the components a = 4-QRq, and an inclination 
angle / = 35°. For simplicity, we assume that the orbit is cir- 



cular and the normal companion is spherical. For a similar 
geometry, Basko (1978) predicted the equivalent width of the 
reflected iron line at inferior conjunction W^^o.s = 7 eV. Basko 
assumed Solar chemical composition of the normal compan- 
ion. However, some X-ray observations indicate that the iron 
abundance may be twice Solar (Kitamoto et al. 1984). Such an 
overabundance of iron would increase the fluorescent line flux 
to W w 14 eV. A narrow iron Ka line of similar amplitude was 
indeed observed by the ASCA SIS (Ebisawa et al. 1996). We 
adopt W0=o.5 = lOeV, which agrees with both theoretical pre- 
dictions and current observations. 

As was discussed in § g, we assume that the Ka light curve 
is measured in the 20 eV energy interval centered at 6.398 keV, 
and the continuum light curve is observed in the 7.1-9.1 keV 
band. The continuum spectrum is assumed to be a power law 
with a photon index F = 1.65 (Liang & Nolan 1984). We also 
allow for the possible presence of a broad iron line originating 
in the vicinity of the black hole, for example, in the outer part 
of the accretion disk. We assume that the equivalent width of 
the broad line is 0.1 keV (see, e.g.. Dove et al. 1998) and that 
the line width is 0.08 keV (FWHM) which corresponds to the 
velocity dispersion of 2000 km s~' . 

For the source variability, we adopt the shot noise model of 
Lochner et al. (1991). This model describes the variability of 
Cyg X-1 as a constant component plus a sequence of exponen- 
tial flares with durations distributed in the 0.01-3.22 s interval. 
The variability level during the 1997 HEAO 1 A-2 observation 
modeled by Lochner et al. was quite typical of the source (35% 
rms). The shot noise model is adequate for simulation of the 
direct emission from the black hole. The light curve of the re- 
flected radiation is the convolution of the direct emission light 
curve with the reflection response to the instantaneous flare. 
This response function is calculated below. 



Black Hole 
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Fig. 1 . — Geometry of the reflection problem. The line of sight is along the 
z-axis. Black hole is in the y-z plane; ip = — ;cos27r0, where is the orbital 
phase. 

Consider the fluorescent emission coming to the observer 
from the normal companion surface element dS (Fig. |l]). The 
distance between the surface element and the black hole is 



I = [R^ + a^-2Ra(cosipsm6cosip+simpcos( 



1/2 



(5) 



The reflected emission reaches the observer with a delay r 

CT = l + asinip-R cos 6 (6) 



See http://constellation.gsfc.nasa.gov; for information 



VIKHLININ 



0.6 



= 0.5 



<t> = 0.25 



= 



{a-2R)/c {a-R)/c a/c 

Time delay 



{a + R)/c 



Fig. 2. — Time response of the fluorescent emission in tlie Cyg X-1 system 
(R/c = 40 s and a/c = 93 s in tliis case). 
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Fig. 3. — The cross-correlation analysis of simulated 1000 ks observations 
of Cyg X-1 at three orbital phases. Correlation functions for = 0.5 and = 
are offset for clarity by 0.003 and —0.003, respectively. 



relative to the direct emission from the black hole. The fraction 
of the X-ray luminosity intercepted by the surface element is 

dL = LdS cosQ il Ait f-, (7) 

where 8, is the incident angle of the X-rays, 

cos 8; = (cos if) smO cosLp + sm if) COS 0)a / 1 - R / d . (8) 

The ratio of intensities of the observed Ka line and the inci- 
dent continuum above the K-edge is a strong function of 8, 
and the fluorescence angle, 6. This ratio, 7(8,, 6*), can be cal- 
culated analytically following the work of Basko. It was also 
derived by George & Fabian (1991) using Monte-Carlo simula- 
tions. George & Fabian present }'(8,, 9) for an incident power- 
law spectrum with F = 1 .7 on a grid of 8, and values. We use 
numerical interpolation of the George & Fabian results in our 
calculations below. 

The light curve of the fluorescent emission from the surface 
element is 

m = Y{Qi,e)5{t-T)dL (9) 

Integration of equation (^ over the normal companion's surface 
provides the time response of the reflected emission. Figure g 
shows the results for the Cyg X- 1 system at three orbital phases. 
A different set of responses is expected for a non-spherical nor- 
mal companion, for example, the one filling its Roche lobe. 
However, the calculations in this case should be analogous to 
the procedure outlined above. 

Our simulations proceed as follows. The continuum light 
curve in the 7.1-9.1 keV band is simulated using the shot noise 
model of Lochner et al. The reflected Koc line light curve is 
the convolution of the 7.1-9.1 keV light curve with the time 
response from Fig g. The Ka line flux is rescaled so that the 
line equivalent width (relative to the F = 1 .65 power law con- 
tinuum) is 10 eV at = 0.5. The 6.388-6.408 keV flux as a 
sum of the Ka line flux and the contribution of the direct emis- 
sion calculated as the 7.1-9.1 keV light curve scaled according 
to the spectral model which consists of the power law and a 



broad 6.4 keV line (see above). Finally, we scale the simulated 
light curves to reproduce the observed Cyg X-1 flux at 7 keV 
(0.05 ph s"' cm"^ keV"', Ebisawa et al. 1996), multiply them 
by the effective area of 1000 cm^, and add Poisson noise. 

The cross-correlation functions derived from the simulated 
lOOOksec observations are presented in Fig 51 The secondary 
peak due to reflected emission is clearly visible at orbital phases 
4) = 0.5 and 0.25. For (f) = Q, the secondary cross-correlation 
peak is weak and not separated from the autocorrelation of the 
direct emission (the average delay is only 14 s for the adopted 
system parameters). 

4. DISCUSSION 

We have shown that the time delay between the direct and 
reflected emission can be measured in a Cyg X-1 like sys- 
tem. Arguably, this measurement requires long observations. 
For a convincing detection, a 250 ksec observation on a tele- 
scope with 1000 cm- effective area around 6 keV is required. 
A longer, 1000 ksec, exposure is needed for a quality measure- 
ment of the cross-correlation function (Fig |]). These exposures 
are longer than the 5.6 day orbital period of Cyg X-1, and there- 
fore the observation needs to be split into several intervals, each 
performed at the same orbital phase. The required exposure is 
shorter for larger area telescopes. Unfortunately, the gain is 
smaller than a factor of (area)'/^ because the internal source 
variability becomes the dominant source of noise. To achieve 
a detection significance similar to that in Fig. || « 300 ksec 
exposure is a minimum even for very large area telescopes. 

In addition to the fluorescence from the surface of the normal 
companion, iron Ka lines can form in the accretion disk or stel- 
lar wind. These emissions are contaminants for the time delay 
measurement. Fortunately, it may be possible to separate them 
in the energy or time domains from the normal component fluo- 
rescence. First, lines originating in the accretion disk or stellar 
wind are likely to be significantly broadened because the re- 
flecting material moves at high speed (a possible detection of 
the broadened iron line in Cyg X- 1 is reported by Done & Zy- 
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cki 1999). Thus, we can exclude a large fraction of the disk and 
possibly wind fluorescent emission by selecting a narrow range 
around the Ka line rest frame energy. Second, iron in both 
wind and disk is likely to be partially ionized because of the 
low matter density and high ionizing flux from the black hole. 
Ka lines for ionization states higher than FeXV are separated 
from the Pel lines by more than 20 eV (House 1969). Third, 
the accretion disk or stellar wind fluorescence should produce 
cross-correlation functions very different from those in Fig. ^ 
Since the outer disk radius is only a small fraction of the black 
hole Roche lobe size (Shapiro & Lightman 1976), the time de- 
lay of the fluorescent disk emission is much shorter than that 
in Fig. ra. In the case of stellar wind fluorescence, time delays 
should be distributed in the broad range At « r/c around the 
average value r = r/c, where r is the characteristic radius of 
the fluorescence region. Therefore, a very broad feature in the 
cross-correlation function is expected due to stellar wind fluo- 
rescence. Such a feature should be easily distinguishable from 
the relatively narrow peak arising from the normal companion 
reflection. To conclude, the only consequence of the presence 
of the iron K-fluorescence from the stellar wind or accretion 
disk is to add noise to the cross-correlation function in the time 
delay range of interest. 

The quality of the data in our simulated observation is suf- 
ficient to determine both the average delay and shape of the 
cross-correlation peak. The average delay provides the orbital 
separation between the components which can be used to refine 
the black hole mass estimates (§ [l]). The width of the correla- 
tion peak is proportional to the radius of the normal compan- 
ion. Its shape depends on the shape of the companion surface 
and the inclination angle (e.g., the (j) = 0.25 curve in Fig. cor- 
responds to all orbital phases in an / = system). The shift 
of the average delay as a function of the orbital phase also de- 
pends on the inclination angle and the companion radius (eq. g. 



Fig. ||). Therefore, a much more detailed modeling than simple 
considerations outlined in eqs. (g)-(Q) might be possible with 
high-quality time delay data. This should result in the compact 
object mass determination without the sin / or Mopt degeneracies 
ofeq(|). 

Basko (1978) pointed out that the orbital modulation of the 
reflected Ka line flux can be used to estimate the orbital param- 
eters, in particular, the inclination angle. Our method is inde- 
pendent of Basko's, because we do not use the line flux. The 
interpretation of the orbital modulation of the Ka line flux can 
be uncertain. For example, the fluorescence in the stellar wind 
reduces the orbital modulation of the line flux; a similar behav- 
ior is expected for low inclination angles (see Basko's Fig. 4). 
On the contrary, our method is not directly affected by either 
stellar wind or accretion disk fluorescence. 

The time delay technique opens the possibility of mass de- 
termination in X-ray systems without known optical counter- 
parts, such as bright black hole candidates in the Galactic cen- 
ter region or LMC. In these systems, the binary period and the 
epoch of zero orbital phase can be determined by periodic vari- 
ations of either the iron line intensity (Basko 1978) or location 
of the cross-correlation peak. Near the orbital phase = 0.25, 
the cross-correlation function is independent of the inclination 
angle and can be used to measure the orbital radius, a, and the 
radius of the normal companion, R, (Fig. 0). The total system 
mass Mopt + Mv is then determined from the Kepler's third law. 
Under the assumption that the normal companion almost fills 
its Roche lobe, which should be the case in X-ray luminous 
systems, the ratio R/a defines Mj/Mopt and thus provides an 
estimate of M^. 
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